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Abstract 

I survey the use of the Haag expansion as a technique to solve quantum field 
theories. After an exposition of the asymptotic condition and the Haag expansion, I 
report the results of applying the Haag expansion to several quantum field theories, 
including galilean-invariant theories, matter at finite temperature (using the BCS 
model of superconductivity as an illustrative example), the Nambu-Jona-Lasinio 
model and the Schwinger model. I conclude with the outlook for further development 
of this method. 
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1. ADVERTISEMENT 

Recently there has been a renewal of interest in analyzing quantum chromody- 
namics starting from the action of the theory and using continuum methods, rather 
than the lattice methods that have been extensively pursued. Light-front methods 
have been emphasized as having the following virtues: 

• Simple vacuum structure 

• Simple boosts 

• Intuitive wave function picture 

• Can be systematically improved 

• Nonperturbative 

Light-front method have some drawbacks: 

• Not explicitly Lorentz invariant, so rotations are complicated 

• Functions are required in counter terms for renormalization 

• The gauge is fixed, so gauge invariance is difficult to check 

• Presence of zero modes 

The amplitudes in the Haag expansion M share the virtues of the light front method, 
and in addition have the following virtues: 

• Have the same number of kinematic variables as Schrodinger amplitudes 

• Obey three-dimensional equations that are explicitly covariant 

• In other words, close to completely on-shell as possible in field theory 

• Composite particles are treated in parallel with elementary particles 

• Can be made crossing symmetric 

This last property holds because, unlike the Tamm-Dancoff expansion, the Haag ex- 
pansion is made in the fields, rather than in the states. The Haag expansion method 
also has some drawbacks: 

• There are more graphs, because many of the lines are on-shell 

• It is not clear in general how to truncate the expansion 

• In confining theories, a replacement must be found for asymptotic fields 
There are some hints how to do this from the Schwinger model. If you like, you 
can use the Haag expansion on the light front and thus combine the virtues and 
drawbacks of the two methods. 
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2. ASYMPTOTIC FIELDS AND THE HAAG EXPANSION 

Asymptotic fields have been part of quantum field theory at least since the 
work of Lehmann, Symanzik and Zimmermannp|; however, there are still miscon- 
ceptions that should be cleared up. The asymptotic fields are free fields, because (at 
least in theories which have neither massless particles nor confinement) the particles 
described by the asymptotic fields separate for large magnitude of the time. The 
free field property of the asymptotic fields does not depend on unphysical "adiabatic 
switching off" of the interactions. The physical assumption is that for t — ► ±00, 
particles either (a) separate widely and thus move freely, since interactions fall off 
exponentially in space or (b) stay close together and thus form a bound state that 
itself moves freely. In this case asymptotic fields must be introduced for the bound 
state. In either case, the exact eigenstates can be labelled by the quantum numbers 
of free particles. The limits for t — > ±00 are the out or in fields that make eigen- 
states at the corresponding limiting times. The asymptotic fields at finite times are 
the limiting fields brought back to finite times according to the free equations of 
motion. The unitary relation between these fields is given by the S-operator, 

S(f) out S ] = cf> in . (1) 

We need asymptotic fields for those bound states that are stable in the approxi- 
mation under consideration; for example, in considering strong interactions, pions 
would be taken to be stable and would receive asymptotic fields. 

The in fields all have free commutation or anticommutation relations and free 
equations of motion and commute or anticommute among themselves. The same is 
true for the out fields. Given the masses and spins of the fields, each set, in or out, is 
a completely known set of fields. Thus each set is a convenient set of building blocks 
for the construction of solutions of quantum field theories. The relation between the 
in and the out fields is nontrivial, given by the S-operator. 

The limits that define the asymptotic fields are subtle. The relations that 
appears in some books, 

<f){x) -> <p out ' in (x), x° -> ±00 (2) 

are ill-defined. The proper limit is a weak operator limit that constructs an asymp- 
totic field of a given mass m from the neighborhood of the mass m part of the 
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relevant (product of) (scalar) Lagrangian field (s) , 

<P° ut > m (x) = lim T ^ ±00 [- f A(x - y; m 2 ) d y o <P(y)d 3 y]. (3) 

JyO=T 

The relation to the mass m part of the Lagrangian field is transparent from the 
momentum space version, 

0<mt, m( k ) S ( k 2 _ m 2) = Um T ^ ±oo ^k )5(k 2 - m 2 ) J dq°(k° + g°)0(g°, k)e- i(9 °- fc ° )r . 

(4) 

For composite particles, one must use a product of the Lagrangian fields of the 
elementary constituents. More about that later. When the field strength renor- 
malization diverges, which is generally the case in relativistic theories, one must 
introduce an averaging over time in the definition of the limit. See|3j for details of 
that. 

To motivate the Haag expansion, recall that we expect that a quantum field 
theory of particles in Hilbert space has three complete and irreducible sets of field 
operators. Here complete means that any state in the Hilbert space can be ap- 
proximated by polynomials in the smeared fields acting on a cyclic vector, usually 
the vacuum state, and irreducible means that any operator that commutes with an 
irreducible set of operators is a multiple of the identity. The first such complete 
and irreducible set is the set of Lagrangian fields, i.e., the fields that appear in the 
Lagrangian and in the action of the theory. The second and third such sets are the 
two sets of asymptotic fields, including fields for bound states, if there are any. Since 
by themselves the set of in fields are completely known, they are standard building 
blocks from which the Lagrangian fields can be constructed. The same, of course, is 
true for the set of out fields. The Haag expansion is just the expression of this idea. 
For a theory with a single (scalar) field and no bound states, the Haag expansion is 

oo - 

A(x) = A m (x) + £ / d 4 x, ■ ■ ■ d A x n f n \x - Xl ,---,x-x n ): A^ n \ Xl ) ■ ■ ■ A^\x n ) :, 

n=2 J 

(5) 

where the double dots indicate normal ordering of the in fields [|1], |J. If there are 
bound states, then in fields for the bound states have to be introduced where- 
ever the conservation laws of the theory allow |J. The physical vacuum is the state 
annihilated by the positive frequency parts of the in fields; only the physical vacuum 
enters and it is a structureless state in this formulation. The equations for the 
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(scalar) in fields are 

(□ + m 2 )A^ n \x) = 0, (6) 

[AW(x),AW(y)U = iA(x-y;m 2 ), (7) 

and the different in fields completely commute or anticommute depending as usual 
on whether the fields are bosons or fermions. If spontaneous symmetry breaking 
occurs, the Haag expansion starts with a c-no. term which is the vacuum matrix 
element of the scalar field]?]]. Haag introduced this expansion in 1955 to discuss 
questions of principle. The present effort is aimed a developing a practical calcula- 
tional method based on his expansion. The /w's ("Haag" amplitudes) are multiple 
retarded commutator functions with all but one leg on- shell J?|]. They automatically 
correspond to connected graphs. They obey the Klein-Gordon equation (for scalar 
in fields) in each argument. Because of this, the terms in the expansion can be 
replaced by 

J d 3 x t f^(- ■■ ,x - x t , ■■■)■■■ f x o A^ n \x t ) ■■■ : (8) 

which illustrates that the Haag amplitudes are both three-dimensional and covariant. 
The asymptotic limit applied to A(x) gives (x) in terms of A^ n \x). The 

relation between the (anti) commutators of these two, 

[A^ t \x),A^ t \y)] T = [A^\x),A^\y)] T (9) 

gives unitarity for all processes. The equal-time commutation relations, 

[A(x,t),A(y,t)]_ = 0,[A(x,t),i(y,t)]_ =^ 3 - 1 5(x,y), (10) 

give generalizations of unitarity. (The choice of A as the canonical conjugate is valid 
for theories without derivative coupling.) The Haag amplitudes for bound states are 
like Schrodinger wave functions for the bound states. There are no relative-time 
coordinates as in the Bethe-Salpeter amplitudes. The Haag amplitudes for elastic 
scattering are the scattering amplitudes with one leg off-shell and the other three 
legs on-shell. The Haag amplitudes for higher processes are related to scattering 
and production amplitudes in a more complicated way. This formalism is as close to 
being completely on-shell as is possible in field theory. Among the good features of 
this approach is the fact that the Haag amplitudes have better ultraviolet behavior 
than the totally off-shell time-ordered amplitudes. This is because the latter contain 
phase-space integrals that grow rapidly for large numbers of particles in intermediate 
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states. The Haag amplitudes for the four-dimensional derivative coupling model, 
which is exactly solvable, illustrates this difference @]. I emphasize that there is 
nothing unorthodox about the Haag expansion; what is surprizing about the Haag 
expansion is that, up to now, it has not been developed into a powerful calculation 
method. That development is the goal of the present work. 

3. NONRELATIVISTIC FIELD THEORY WITH BOUND STATES 
3a. The model 

Consider a model with two spinless nonrelativistic Fermi fields, A(x) and B(x), 
with the Hamiltonian 

H = m A [ d 3 x AUx) A(x) -\ — [ d 3 x V x A f (x) ■ V^A(x) + (B terms) 

J 2m A J 

+ / d 3 xd 3 yB\y)A\x)V AB (\x-y\)A(x)B{y); (11) 

Jx°=y° 

for simplicity I assumed an AB interaction, but no AA or BB interaction. The 
equation of motion for A(x) is 

id a oA(x) = (m A --?—V 2 x )A(x)+ I d 3 yB\y)V AB {\x-y\)B(y)A(x). (12) 

2m A Jx°=y 

The asymptotic (in or out) fields for (possibly composite) particles are character- 
ized by their rest energy E, kinetic mass m and spin J. I suppress the spin in what 
follows. The kinetic mass is the mass that enters in the kinetic energy, p 2 /2m; for 
composite particles, as discussed below, the kinetic mass is the sum of the kinetic 
masses of the constituents, without the binding energy, because of the Bargmann 
mass superselection rule described in Sec. 3b. The asymptotic fields obey the fol- 
lowing free field equation and anticommutation or commutation relations: 

id x oC in(out \x) = (E- —V 2 )C in{out \x) (13) 

2m 

[c in(out) (X; t); c \in{out) (y? ^ )]± = p(x _ y> f _ ff . ^ m)> (M) 

D(x, t ;E,m) = ^-f d 4 kS(k° -E- ^-) e -^ t+ ^\ (15) 
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Note that 



V{x,0;E,m) = <5(x), VE,m. 



(16) 



Using translation invariance, on can show that the Haag expansion of the 
interacting field A(x) in terms of in fields takes the following form in position space 
(with an analogous expansion for the B field) |§ 

A(x A ) = A in (x A ) + Yl[ d 3 x B d 3 Xif B]i (x A - x B ; x A - x i )B^ in (x B )(ABi) in (x i ) 

i 

+ J d 3 x B d 3 y A d 3 y B f B . AB {x A - x B ; x A - y A , x A - y B )B ]in {x B )A in {y A )B in {y B ) 

+ •••, (17) 

where, since both the asymptotic fields and the Haag amplitudes obey free equations, 
the integrals are independent of the times x° B , y° A and y B because of the translation 
invariance of the Schrodinger scalar products. Label the Haag amplitude that is the 
coefficient of a product of (asymptotic) creation and annihilation operators by the 
labels of the operators. Labels to the left of the semicolon are creation operators, 
to the right are annihilation operators; the two-body {AB) bound state in level i is 
labeled by i. 

Some calculations are simpler in momentum space, therefore define 

A(x) = J d 4 ke- ik x A{k), (18) 

V AB (\x-y\) = j^y 3 J d 3 qe^-^V AB (ci). (19) 

From now on, drop the tildes on A and V. Transforming the equation of motion to 
momentum space yields 

(k° A -m A - )A(k A ) = / d 4 k B d 4 p B d 4 p A S(k A + k B -p B - p A ) 

2m A J 

xB\k B )V AB (p B - k B )B(p B )A(p A ). (20) 

Define 

C m (x) = (2vr)- 3/2 J d A k5(k° -E c - k 2 /2m c )e- ifc ° t+ik x c in (k), (21) 
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[c m (k), C t-(l)] + = 5(k-l) (22) 

fsMv) = 7^8 / ^^^^e^^^^^^^^-^^^^^-^B^k!,^) 

(23) 

and similar definitions for other Fourier transforms chosen so that powers of 2tt are 
absent from most of the momentum-space formulas. The result is 

A(k A ) = (2n)-^:k A :6(k A -m A - k ' ] > 
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+ / d 3 k B d 3 ki 5(k° A + — -m A + 6i- —^)5(k A + k B - ki)f B;i (k B ; ki) : fc^ : 

l 2 2 2 

+ / d 3 k B d 3 PB d 3 PA 5(k° A + T ^-~m A -^- ^)5(k A + k B - PA - p B ) 
J Zm B Zm A Zm B 



x fB-,AB(k B ; pa, Pb) : k B p A p B : H , (24) 

where I define : p A := a m (p), with normalization as given in Eq. ([22]) , etc. Note that 
I am expanding in terms of in fields; there are analogous expansions in terms of 
out fields. In the next section I derive the constraints on the /'s that follow from 
Galilean invariance. 



3b. Galilean Invariance 



Bargmann||12|| showed that the unitary projective representations (i.e., repre- 
sentations up to a factor) of the Galilean group that occur in the quantum mechanics 
of nonrelativistic particles cannot be reduced to vector (i.e., true) representations. 
This contrasts with the corresponding situation for the Poincare and Lorentz groups 
(and indeed most other physically interesting groups), where the representations can 
be reduced to true representations. The explicit mass parameter in the phases leads 
to the Bargmann superselection rule that the sum of the masses (that appear in the 
kinetic terms) must be conserved in every process. Nonetheless, bound states can 
be formed and particles can be created and annihilated, provided the Bargmann 
superselection rule is obeyed. 

Note that, for example, a bound state of particles of masses m A and m B with 
binding energy e has energy E = m A + m B — e + k 2 /2m AB , rather than E = 
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m A + tub — e + k 2 /2(m^£ — e) as one might expect from the nonrelativistic limit 
of a relativistic bound state with rest energy m A + m>B — e. (I use the abbreviation 
tyiab = rriA + itib-) Another manifestation of this effect is that for this same bound 
state the momentum transforms under Galilean boosts as k — > k + m^v, rather 
than as k — > k + (tjiab — e) v - 

If the projective representation has the form 

U(G 2 )U(G 1 ) = u(G 2 , G l )U(G 2 G 1 ) (25) 

then another projective representation is equivalent to this if the other representation 
has the factor system u/(G 2 ,G'i) = [(f)(G2)4>(Gi) / (f>(G2Gi)]uj(G2, Gi) , where <p has 
modulus one. This arbitrariness allows simplification of some formulas. 

Bargmann gives as the Galilean transformation of a nonrelativistic scalar wave 
function, 

(T(G»(x) = e- id ^ G ' x) i){G~ l x), (26) 

where x = (x, t), the Galilean transformation is Gx = (-Rx + vt + a, t + b), G = 
(a, b, R, v), where a and b are space and time translations, R is a rotation and v is 
a boost, and 6{G,x) = m(|v 2 t — v • x). To infer the corresponding transformation 
for a nonrelativistic scalar field, I require 

U(G)A(ip)U\G) = A^g), Vto(aO = (T{G)^)(x) = e-^V(G~H (27) 

A($) = J A(x)^(x)d 4 x. (28) 

Then 

U{G)A{x)U\G) = e' idA(G ' Gx) A(Gx), 9 A (G,Gx) = m A [^v 2 (t+b)-v(Rx+vt+a)}. 

2 (29) 

If the field has spin s, then A on the left hand side is replaced by Aj and A on the 
right hand side is replaced by J2j AjDjf(G), where is a representation of SU (2), 
which is the little group in this case. The corresponding transformation holds for B 
with rriB replacing m A . Asymptotic fields transform the same way. The implications 
of the transformation law for the Haag amplitudes is found by transforming the 
interacting field in two ways: (1) insert the right-hand side of Eq.(p9|) in the Haag 
expansion, or (2) transform each of the asymptotic fields and then change variables 
to get the transformation into the amplitudes. The two amplitudes fs-i and fs-AB 
obey 

fB-A G ( x A - xb)\ G(x a - Xi)) = 
9 



e ie A (G,Gx A )+ie B (G,G XB )-ie AB (G,Gx i ) fB .( XA _ Xb . Xa _ (3Q) 
fB-,AB(G(x A - x B )] G(x a - y A ), (G(x A - Vb)) = 

Note that 9 AB is independent of the bound state i because of the Bargmann mass 
superselection rule. The combination of phases in the first of these is 

6 A (G, Gx A ) + 6 B (G, Gx B ) - Oab(G, Gxi) = 

- ^v 2 (m A x° A + m B x° B - rrijisx^) - v • R(m A x A + m B x B - m AB Xi). (32) 

The transformation law is not satisfied by having a delta function in the space and 
time coordinates identifying the coordinate Xi with the center-of-mass of particles 
A and B, although at equal times such a delta function does occur for the space 
coordinates. The way in which the transformation laws are satisfied is best seen in 
momentum space, to which I now turn. 

The corresponding transformations in momentum space are 

{V{G)4>){k) = e-^^^G-'k), (33) 



Q(G, k) = (k- mv) • a — (k° — ^mv 2 )&, (34) 

where k = (k, E), Gk = (Rk + mv, E + v • Rk + |mv 2 ), and G _1 k = (ir x (k - 
mv), E — k • v + |mv 2 ). The momentum space transformation law for the field is 
induced in parallel with the derivation of the position space law. The result is 

W{G)A(k)W\G) = e- inA{G '~ Gk) A(Gk), (35) 

where £l A (G, —Gk) = (E + v- Rk)b — Rk ■ a. In the transformation law for the Haag 
amplitudes, all the phase factors cancel and the result for-say-the second term in 
the Haag expansion is what one would expect naively, 

f B -i(k B ; kj) = f B -i(R(k B - m B v); Rfc - m AB v)). (36) 

Thus I can choose the v = ki/m AB so that the bound-state momentum vanishes 
and eliminate the second argument of f B -i, 

f B;l (k B ; k,) = / Bii (k B - ^k u 0) = F B -i(k B - -^-ki). (37) 

m AB m AB 
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For the spinless case, F S;i (k) = F B .i(Rk). All these results are exact, valid in any 
Galilean frame. The extension to fields with spin is straightforward. It is worth not- 
ing that the Poincare transformation law in a relativistic theory is simpler than the 
Galilean transformation law because the Bargmann phase is absent for the Poincare 
group. 

Taking account of Galilean invariance, one finds that the position-space Haag 
amplitude is 

y) = (27T)" 3 / d 3 kd%exp[t(m B + -L(k + -^-k,) V - i(k + ^-k,) ■ x] 

J 2m B m AB m AB 

k 2 

xexp[-i(m AB - q + — —)y° + ik* • y]/s ; i(k; 0). (38) 

2m AB 

The integral over kj can be done, but the result is complicated and not useful, except 
when all times are equal, in which case the result is both simple and useful, 
t t \ r/ m A x A + m B x B 

Jb-A^a ~ x B ; x A - Xj) = d(Xj )F B . ti (*A - x B ), (39) 



F B;l (x) = | rf 3 A;e- jk - x /B,(k;0). (40) 



Using the constraints due to Galilean invariance, the Haag expansion in rr-space at 
equal times takes the form 

,m A x. A + m B xg, 
m AB 



A(x A ) = A m (x A ) + W F B ,( XA - XB )B^ B )(AB i y n ( WA + mB * B )d 3 x B 
. J in a r 



+ [ dWrF B . tAB (r'; r)B*»(x A - r>)A™(x A + m *l r T) )B^ A - mAT + m * r ) 
J m AB m AB 

+ •••• (41) 



In momentum space, the expansion is 



1 k 2 

A ( k A) = / n/2 : ^ : 5 ( k A -m A - ' 1 



( 27r )3/2 • • ^ — 2mA J 

frfi w/o k ^ j_ (k A + k B ) 2 ~ m A k B -m B k A f 

+ / d- i k B 5{k A + - m^ + ei )F B;i ( ) : k [ B (k A + k B )i : 

J 2m R 2m a r Wjr 



I 2 2 2 

+ / d 3 k B d 3 p B d 3 p A 5(k A + —^--m A -^A_- ^)5(k A + k B -p A - p B ) 



2m B 2toa 2m B ' 
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~ m A k B -m B k A m A p B - m B p A f 

x F B . AB ( ; ) : k B p A p B : H . (42) 

m AB m AB 

3c. Two-Body Bound State 

To derive the equation for the two-body bound state, insert the Haag expan- 
sion Eq.([T7]) in the equation of motion Eq.([TJ), renormal order, and equate the 
coefficients of the terms with the operators B^ in (ABi) m . After commuting or anti- 
commuting with the relevant in fields, the result is 

d 1 

^daP ~ mA+ 2~m _V ^ A ~ V ^ A ~ x B))fB-,i( x A ~ x B ; x A - x{) = (43) 

It is convenient to eliminate the time derivative by using d/dx° A = —d/dx° B — d/dx®, 
the independence of the Schrodinger scalar product on the time and the free equa- 
tions satisfied by the in fields to find free equations for the x° B and x° dependences 
of f B -i. The results are 

" mB + ^f« )fB * = °' (44) 

(! S? +m '«- f '-2^ V - )/fl »= a (45) 
The equation without time derivatives is 

[-i v - - sb v « + v ^ ~ Xb)1/ «" = (e - " 5^^" (46) 

Now using Eq . (|39|) the usual Schrodinger equation for F B]i results, 

[~^l AB + V(v AB )]F B , = -e^ 1 = -L + -L, (47) 
z/i fx m A m B 

where the reduced mass enters. This establishes that F B -i is the Schrodinger wave 
function of the bound state. Note that the bound-state amplitude is given exactly 
in any reference frame in terms of the amplitude in the rest frame of the bound 
state. (The corresponding statement also holds for other amplitudes, as well as for 
relativistic theories.) 
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3d. Two-Body Scattering 



Two-body scattering is described in position space at equal times by the am- 
plitude 

/B;AB(xA-x B ,0;x A -y A ,0,x B -y s ,0) = F b . ab (x a -xl b ; y A -y B )5(R'-R), (48) 

FbM^y) = (27r)- 3 / 2 ! d 3 k'd 3 kf B , AB (k'-,-k,k)exp[t(-k'i^ A -^ B )+kiy A -y B )], 

(49) 

, m A x A + m B x B m A y A + m B y B 

XV — , rv — . 

m AB m AB 
I prefer to discuss two-body scattering in momentum space, using the amplitude 

fB-,AB(k B ; p A ,p B ) which is the coefficient of the term : k B p A p B : in the Haag expan- 
sion of A(k). The procedure for finding the equation for f B;AB is analogous to that 
for the two-body bound state amplitude. One finds 

( P A -(pA + PB-k B ) 2 p|~k| - 

{ ^ 1 o jB;AB{k Bl p A , p B - 

2m A 2m B 

v AB (\k B - p B \) + J d 3 k' B v AB (\k B -k' B \)f B . )AB (K;PA,PB). (50) 

Galilean invariance relates f B - AB at arbitrary momenta to itself in the center-of- 

mass, 

fB;AB(k B ; pa, Pb) = ]b;Ab{R{^b - m B v); R(p A - m A v), R(p B - m B v)). (51) 

By choosing v = (p^ + p B )/m AB , I can replace f B - AB by a function of one fewer 
variable, 

Ib\AB (k B ] p A , p B ) = F B . AB (k; p), (52) 

where here and below, k = {m A k B — m B k A )/m AB , p = {m A p B — m B p A ) /m AB and 
I used conservation of momentum to introduce k^. The momenta p and k are the 
center-of-mass momenta of particle B in the initial and the final state, respectively. 
The elastic scattering equation becomes 

^(p 2 - k 2 )F B , AB {k; P) = V{\k - p|) + J d 3 k'V(k - k')F B;j4B (k'; p), (53) 

The solution is the Born series, 
F B . AB (k; p) = 

G R (k; p)V(\k - p|) + G R (k; p) J d 3 k'V(\k - k'\)G R (k'; p)V(\k' - p|) + • • • , 
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where G R (k; p) = [(p 2 - k 2 )/2u - ze]" 1 . 

The amplitude F B . AB is closely related to the T-matrix element for AB scat- 
tering. The S'-matrix element is 

S{k A , k B ;p A ,p B ) = out (k B , k A \p A ,p B ) in = (0| : k^k^p^ : |0), (54) 

where I remind the reader that : k A :, etc., stands for the in field. In order to 
eliminate the out fields in terms of the in fields, use the definitions, 

A mi ~ out \x)= lim f d 3 yV{x-y-m A ,m A )A{y), (55) 

T^ — oo(oo) JyO =T 

where T> was defined in Eq.([T5|). The nonrelativistic analog of the reduction formula 
follows from calculating / d^yd / d y oT>(x — y;m A ,m A )A(y) in two ways: performing 
the integral and carrying out the derivative. The result [|J is 

A° ut {x)-A m {x) = J d 4 yV(x-y;m A ,m A )(d y o+tm A -^-V 2 y )A(y). (56) 

Fourier transforming this one gets, after removing a factor of S(k° — m A — k 2 /2m A ), 

1 k 2 

{a out (k) - a in (k)) = -2ni(k° - m A )A(k). (57) 



The right-hand-side is non-vanishing (and there is scattering) only when A(k) has a 
pole at {k°-m A -k 2 /2m A ) = 0. Since a f OU *(A;)|0) = a f in {k)\0) for stable particles, 
the only out operator in the S'-matrix element (0| : k ^ 1 :: k ^ :: p A :: p B n : |0) 
that must be eliminated using Eq.(]57f) is : k™ 1 :. The result is 



S(k A ,k B ;p A , PB ) = 5(k A - PA )5(k B - PB ) -2m5(^ + kjj P:i PL ' 



2771^4 2m B 2m A 2m B 

x ^ + kB - PA - PB)l 3. + ^.-A.- go^fc p), (58) 

where again k and p are defined below Eq.(^). Thus the reduced T-matrix for 
elastic scattering on the momentum shell is 

ks; >*• PB) = ^ + - ^ - H^-* < 59) 
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I emphasize that because the Haag amplitude is the scattering amplitude with one 
leg off shell, it contains the information necessary for calculations in the three-body 
sector. This contrasts with the on-shell scattering amplitude, which does not suffice 
for such calculations. 

3e. Anticommutation Relations 

In this section I show that the canonical (equal time) anticommutation re- 
lations of the Lagrangian fields imply general relations among Haag amplitudes, 
independent of the equations of motion of the specific theory. For example, the 
vanishing of the canonical anticommutator [A, B]+ at equal times, considered for 
the coefficient of the bound state in field for state i, gives 

F A;i {y - x) = F S;i (x - y) = F<(x - y) (60) 

where I took (ABi) m (R) = — (BAi) m (R) because of the Fermi statistics of A and 
B. Thus the apparent asymmetry in the treatment of the constituents of the bound 
state, due to the fact that the Haag amplitude that serves as the two-body wave 
function of the (AB) bound state in the Haag expansion of the A field has the A 
particle off-shell and the B particle on-shell, while these roles are interchanged for 
the amplitude for the same bound state in the Haag expansion of the B field, is 
not a real asymmetry. These two amplitudes determine each other uniquely. The 
analogous result for the off-shell elastic scattering amplitudes is 

F B ;Ab{* - y; r) = F A]BA {y - x; -r) = F AB {x - y; r). (61) 

Again the two apparently different off-shell amplitudes uniquely determine each 
other. 

The consequence for elastic scattering is 
t(k A , k B ; p A , p B ) - (t(p A , p B ; k A , k B ))* = 

(2vr) 5 / 2 / d 3 q A d 3 q B S(^- + J^L - -SL - ^)5(k A + k B - q A - q B ) 
J 2m A 2m B 2m A 2m B 

xt(k A: k B ; q A: q B )(t(p A , p B ; q^, q B ))*, (62) 



where k and p are as defined below Eq . (|52|) . This is elastic unitarity. 
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The canonical anticommutator [A, A^] + at equal times leads to a generalization 
of unitarity, 

j—^(F B . AB (k; p) + F%. AB (p; k)) 
= ^F B;J (k)F* , t (p) + J A^;A B (k;q)F* AS (p;q), (63) 

i 

where again k and p are as defined below Eq.(^) and I have used momentum 
conservation, k^ + k B = Pa + Pb- By taking the appropriate limit, I recover the 
elastic unitarity relation, Eq. (|62"D . On taking into account the relations between the 
Haag amplitudes in the expansions of A and of B, one find that these are all the 
independent two-body relations obtained from the anticommutation relations. 

There are also quadratic relations between the amplitudes for the (ABi) and 
(ABj) bound states and the amplitudes for the breakup of these bound states due 
to scattering with the A or B particle. Since this involves a higher sector, I do not 
give this relation here. 



3f. Construction of the asymptotic field for the bound state 



In this section I show how to construct the asymptotic field for the bound 
state from a product of Lagrangian fields. My suggestion differs from that proposed 
by Nishijima|T(J and by Zimmermann[ll]. The procedure is to multiply the appro- 



priate Lagrangian fields at separated space points, integrate with the bound-state 
amplitude in the relative coordinate, and take the asymptotic limit. If the in field 
expansions of the Lagrangian fields are inserted and the resulting expression normal 
ordered, then the t — > — oo limit gives the in field bound state operator and the 
t — > oo limit gives a reduction formula for the out field bound state operator. The 
result is 



{ABif< out \x) = lim / d 3 yV(x~y;m AB -e l ,m AB )F 

r^ — 00(00) JyO =T 



w) 



x-[B(y-—w),A(y+—w)]d 3 w. (64) 
2 m AB m AB 

A straightforward calculation shows this limit is (ABiy n (x) for r — > — 00 and the 
leading term for r — > 00 is (ABi) out (x). Both results are what we expect. The 
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later terms in the Haag expansion for (ABi) out (x) are in a higher sector that I don't 
discuss here. 

I derived many results of the nonrelativistic quantum mechanics of two-particle 
systems in a unified way with particular attention to Galilean invariance, taking into 
account the fact that the representations of the Galilean group in quantum mechan- 
ics are necessarily representations up to a factor, rather than vector representations. 
The Haag amplitude for the simplest term with the two-body bound-state opera- 
tor is precisely the Schrodinger wave function of the two-body bound state. The 
amplitude for the term with three in fields is the scattering amplitude with one leg 
off-shell. These interpretations carry over to explicitly covariant relativistic theories, 
where the corresponding Haag amplitude is defined on three-dimensional manifolds, 
but is covariant. Of course in the relativistic bound state that is mainly a 

two-body state also will have amplitudes to be composed of more than two particles. 

4. The NAMBU-JONA-LASINIO MODEL 

The Haag expansion is effective in treating the Nambu-Jona-Lasinio model 
in one-loop approximation. In particular, the Haag expansion sums crossed as well 
as direct graphs [|T3|, in contrast to the usual methods which sum only direct 
graphs. Further, using the Haag expansion one deals directly with the bound-state 
wavefunction (or amplitude); one does not have to extract the bound-state amplitude 
as the residue of a pole in the scattering amplitude. 

The Lagrangian of the model without isospin is 

£ = iffij, - - go [(^ 7jU ^) 7 ^ - (^7 M 7 5 ^)7 M 75# ( 65 ) 

To take account of operator ordering, symmetrize or antisymmetrize the operator 
products. After going to momentum space, the equation of motion is 

^(<?) = -~go J d i p 1 d 4 p 2 d*p 3 5(q-pi-p 2 -p3) (66) 
{[[^(pi),V(P2)]-,V'(P3)]+ - Uipi),l^ip2)Ul^{pz)] + } (67) 

The lowest approximation, to take i/j(p) = ip m (p)5(p 2 — m 2 ), leads to the gap equa- 
tion, 

/i = ^/ rfW " /i2) - (68) 



Surely, this is a simple derivation of this result! Since this model is nonrenormal- 
izable, one must cut off the integral. This can be done covariantly, if desired. For 
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< ir 2 /g A 2 < 1, where A is the cutoff, there are three solutions: /i = and 
= ±m. The first is the unbroken symmetry solution; the last two are equivalent 
broken symmetry solutions. To decide which solution is the stable one, calculate the 
matrix element of the Hamiltonian in the corresponding vacuum state. The result 
is that the broken symmetry solutions have lower vacuum matrix elements and are 
thus the correct solutions. 

To find the bound states, consider the term in the Haag expansion with the 
product : if) in B in :, where B m is the bound state in field. The coefficient of this term 
serves as the bound-state wavefunction; indeed as shown above in the nonrelativistic 
case it is precisely the Schrodinger wavefunction. Inserting the two terms of the 
Haag expansion that have been introduced into the equation of motion and re- 
normal ordering and keeping the coefficients of the bound-state term leads to a 
linear integral equation for the bound-state wavefunction. Because the interaction is 
a contact interaction, this equation can be solved exactly. For the J PC = 0~ + state, 
the analog of the pion, the mass is zero, as expected from the Nambu-Goldstone 
theorem. For other states, the results agree generally with previous calculations, 
except in some cases the limits on the masses differ, perhaps because the present 
calculation includes crossed graphs. 

5. FINITE TEMPERATURE FIELD THEORY 



5a. Sketch of thermo field theory using the Haag expansion 



I illustrate the application of the Haag expansion to the solution of second 
quantized field theories at finite temperature using the BCS model of superconduc- 
tivity. In order to have a state annihilated by the annihilation operators, which is 
necessary for normal ordering, I use the thermo field theory formalism of Umezawa 
and collaborators^, |TB[. This formalism uses a doubled set of operators to account 
for finite temperature. The annihilation and creation operators are subjected to 
two Bogoliubov transformations: one comes from the dynamics of the electron pair 
interaction; the other from the thermo field formalism which takes account of the 
finite temperature. In lowest approximation, the method leads to the usual gap 
equation. The asymptotic fields whose annihilation parts annihilate the vacuum at 
zero temperature no longer annihilate the state which is a thermal mixture at finite 
temperature T. Indeed, no set of annihilation operators annihilates the mixed state 
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at finite T. In order to obtain a state which is annihilated by the annihilation parts 
of a suitable set of asymptotic fields, the Hilbert space of states must be enlarged 
to include hole states in the thermal equilibrium state at a given temperature and 
the set of operators must include operators which annihilate and create holes in 
addition to the operators which annihilate and create particles. Thermo field theory 
does this, for example for a Hilbert space with a discrete energy eigenstate basis 
{|n)}, by replacing the Hilbert space of the system under consideration with the 
tensor product Hilbert space with basis {\n) <g> \n)}. Correspondingly the set of op- 
erators with respect to which the vacuum |0) ® |0) is cyclic is doubled and includes 
the particle {a& ® 1, a' k <S> 1} and hole {1 <g> a^, 1 <S> %} operators. (To simplify the 
notation, I will drop the ®1 and 1® factors.) A Bogoliubov transformation relates 
the original annihilation and creation operators for the particles together with the 
"tilde" operators which describe the holes to another doubled set of operators whose 
annihilation parts annihilate the pure state (called the "thermal vacuum") in the 
enlarged Hilbert space which represents the thermal mixture in the usual theory. 

Let the density operator be 

p = Z{fi)- x e- m = Z03)- 1 J2 e~" En \n) (n\, (69) 

Z(P) = tre-P H = Y,e~ f>En - (70) 

Here H can either be the Hamiltonian H or H — pN, where p is the chemical 
potential and N is the number operator. Now consider an enlarged Hilbert space in 
which the tensor product basis {|n) <8> \n)} replaces the basis {\n)} of the original 
Hilbert space. Let the thermal vacuum be 

\0(P)) = Z{P)-^ 2 ^2e-^ 2 \n) ® \n). (71) 

Let the doubled set of operators, c, c' , c and ct be operators for the normal modes 
of the total system. The thermal vacuum at finite temperature is the state which 
satisfies 

c kl \O((3)) = 0, hiW)) = 0. (72) 

The existence of a state \0((3)) which is annihilated by the annihilation operators 
is essential to defining a normal-ordered operator product. The average of an ob- 
servable A in the thermal mixture at inverse temperature (5 is given by the matrix 
element of the corresponding operator A <g> 1 in the thermal vacuum \0((3)), 

tr{ P A) = (0{(3)\A®l\0{(3)). (73) 
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For the Fermi case of interest for the electrons in superconductivity, the Bogoliubov 
transformation between the operators, b and b\ for the normal modes of the electrons 
in the system and the operators, b and &t, for the normal modes of the holes on the 
one hand and the doubled set of operators, c, c', c and d" which are the normal 
modes of the total system on the other hand is 




\/l - n k Jnk~ 



kl 



n k y/l - n k J \ c k i 



V 1 - n k -\frik \ I c fc2 
b{ 2 ] \ V™k V±-nk I I 4 2 



(74) 



(75) 



The requirement that 

(0(P)\blb kl \Om = (0(P)\bl 2 b k2 \0((3)) =n k = l/(eP B - + 1) (76) 



fixes the coefficients in the Bogoliubov transformation. As indicated in Eq. (76), I 
assume that n k is independent of spin polarization. 

The electron-electron interaction which leads to superconductivity leads to 
another Bogoliubov transformation which has the form 



4i 

a k[ 



a k] 




(77) 



and 



o-fct 1 I u k 



a k 




(7f 



The solution of the operator equations of motion for the a and a operators determine 
the coefficients u and v in Eq. flTTD; the operators b and b are then expressed in 
terms of the c and c set. The final result gives the a and a operators in terms of the 
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c and c operators, 



( 4 T N 

a -k[ 

V G-*4 / 



-v 



l^Jl - n k u* k y/l - n k -v* k ^Jnk~ 



-v k ^n k 



-u ky /n k 
v* k \/nk 



v ky /n k 



u k \/l - n k v k y/l - n k 



-VfcVl - n k u* kV /l - n k J 



1 c ki ' 

Cfc2 
V 4 / 



(79) 



5b. Solution of the finite temperature equations of motion 

Assume there is an attractive interaction between pairs of electrons with op- 



posite spin and momenta which acts only on electrons near the Fermi surface |T7 
The total Hamiltonian is 

H = H — H, (80) 
H = T, ks (k 2 /2m)a\ s a ks + (1/4)£ M [<4 T , ai kl ]V ke [a-ei, a £T ], (81) 
where the potential acts only near the Fermi level, 



Vq, \h 2 p 2 /2m — fi\ < Ae, p=k or I 



Vi 



ki 



0. 



otherwise. 



(82) 



The tilde Hamiltonian, H, is the same as H, except that a and a+ replace a and a'. 

It is convenient to measure energy relative to the Fermi surface and to use 
H — fiN instead of H as the energy. This has the effect of changing the coefficient 
of the kinetic term in the Hamiltonian from k 2 /2m to e k = k 2 /2m — \i. I retain 
the symbol H for this new "Hamiltonian" and the symbol H for the new "total 
Hamiltonian." 

The time dependence of any operator is 0(t) = e %Ht Oe~ lHt . Thus the equation 
of motion is 

-id t O(t) = [H,0(t)). (83) 
The equations of motion for creation and annihilation operators are 



- id t a ] k] = e k a) k] + (l/2)E^[a] T , aL^]V^ fc a_fe|, 

- id t a-ki = _ 4a-fci + (l/2)S^4 T Vj w [a_^, a^}, 

- id t a k ^ = £ k a k] + (l/2)E e a ] _ H V M [a_ ei ,a ei }, 



J4) 

J5) 
i6) 
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- id t a ] __ kl = -e k aL H + (l/2)E < [a{ r , al n ]V ek a kV (87) 

Use a Haag expansion in the thermal in (or out) fields to get an approximate solution 
of the operator equations of motion. In the lowest approximation, use Eq.([79|), insert 
it in the equations of motion, re-normal order and keep only the linear terms in the 
c operators to get equations for the unknown coefficients u k and v k . The result has 
the form 

^(F fe + ^G fe/ )^ = 0, (88) 
i i 

where Cj = (cl 1 ,Ck2,Cki,c\. 2 ) (as a column vector). Since the c and c operators are 
linearly dependent, each of the 16 equations, 

(F k + Gkih = 0, i, j = 1 to 4, (89) 

must hold. Only two of these equations, 

(E k - 4 - V kk w k )u k + T,(VuytV k = 0, (90) 
(E k + £ k + V kk w k )v k + Y^fYuyiUk = 0, (91) 

where 

w k =vl(l -n k ) +u 2 k n k , yi = uivi(l-2ni), (92) 

are linearly independent. It is convenient to write these equations in matrix form 
and to make the equations appear simpler by introducing new symbols 

efe = e k + V kk w k , A k = ^2 v kWh (93) 



Then 

= 0. (94) 



E k -e k A fc \ / u k 
A k E k + e k I \ v k 



There is a solution only if the determinant of the matrix vanishes, i.e., if 

El = e 2 k + Al (95) 



The solution, which is implicit because e and A depend on u k and v k via Eq. (|92| 
and (0), for u k and v k is 



u k = yJ(E k + e k )/2E k , v k = -A k /pE k (E k + e k ). (96) 
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It is easy to check that this solution satisfies the constraint u 2 k + v\ = 1. When 
Eg- (Fop is inserted in the definition of A, the celebrated gap equation results, 

A k = '£V kl {A l /2E l ){l-2n l ). (97) 



With the choices, 



. . , -V , \e k \ < Ae, p = k or I 

Vki — \ 

0, otherwise, 



A k = { A(T) ' |6fe| - Ae ' (99) 
0, otherwise, 

using Eq. ([f6|) and making the usual replacement of the sum by an integral, the gap 
equation becomes 

V N(0) f £ de l tanhipEi^/E^l, (100) 
Jo 



where E k = + A 2 (T). This is the usual gap equation for this model and the 
usual results for A(0)/fcT c = 1.764, N(0)V Zn(1.13Ae/JfcT c ) = 1, etc., where T c is 
the transition temperature, follow. 

6. THE SCHWINGER MODEL 

The Schwinger model |TR [H| is massless two-dimensional quantum electrody- 
namics, an exactly soluble model. In the Lorentz gauge,| 

£ = $i$kl> - - \{d ■ A) 2 , (101) 

where 

D ft = &*-eA ,i , F» v = d»A v -d v A». (102) 



Lowenstein and Swieca|20| found an operator ansatz that yields the matrix 
elements computed by Schwinger. Their ansatz solution has the following form: 

A»(x) = - JJ(e^^S + d*7i), (103) 



2 This form of the Lagrangian ignores wavefunction renormalization of the spinor field. Taking 
account of this wavefunction renormalization, the spinor term in the Lagrangian is ip u ipip u — 
(i/j u 'i'P>Pu)o = Zipipip, where tp u is the unrenormalized spinor field and Z is the spinor wavefunction 
renormalization. We suppress the factor Z below. 
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ij)(x) =: exp[iv^F(7 5 S(x) - T](x))} : ^ (0) (x), (104) 

where rj is a free neutral massless field with negative metric corresponding to the 
gauge degrees of freedom satisfying [rj(x),fj(y)]ET = —io~(x l — y l ) , while S is a 
free neutral massive field with positive metric representing the physical degrees of 
freedom satisfying [T,(x),T,(y)]ET = i6(x 1 — y 1 ) and ip^ is a solution of the free 
massless Dirac equation. (ET stands for equal time.) This solution displays the 
main property of the Schwinger model: the only physical state is a free particle of 
mass e/ s/tt, but the solution does not obey the canonical commutation relations for 
Ap. The Haag expansion provides a solution that does obey these relations [p2fl. 

Use the Lagrangian Eq. (|101| ), but drop surface terms so that the gauge field 
part becomes 

-\{d ll A v ){d> i A v ). (105) 

Since the Lorentz group (without inversions) is abelian in 1+1, all irreducible repre- 
sentations are one-dimensional; thus the vector and spinor fields as to in the model 
are composed of one-dimensional irreducibles arbitrarily pasted together. Express 
the Lagrangian in terms of the irreducible fields in the basis in which 

A° = ±(A + + A-), A 1 = ±(4+ -A"), tf = (Vi,lfcO, (106) 

with 



7 5 = I I " )• (107) 



In terms of the irreducible fields, 

£ = i/;l{2id + - eA")Vi + ^l{2id~ - eA+)^ 2 + )^{d 1 A+d x A- - d°A + d°A-). (108) 

using lightcone coordinates, x + = x + x , x~ = x° — x . The corresponding 
derivatives are = \{-^ ± ^fr). Define these so that -^x ± = 1. Note that 
although the fields A ± are lightcone fields, I am not using lightcone quantization, but 
rather am using equal-time canonical quantization. The naive operator equations of 
motion are 

UA + - 2eV>lV>i = 0, (109) 

UA~ — 2eiplip 2 — 0, (110) 

(2id + - eA~)i) X = 0, (111) 
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(2id~ - eA + )^ 2 = 0. (112) 

As Schwinger pointed out in his original paper, the spinor bilinear products require 
a line integral of the "vector" potential between the ip' and the ijj in order to ensure 
gauge invariance; this is done explicitly below using point-splitting. For example, 
Vwi is replaced by 

1 



lim 

e^O 2 



4(x + e)e~ ie ^ '^^"^(x) + cc. . (113) 
The canonical momenta are 



ir A+ = ~df ) A-, n A - = ~d°A + , 7r^,=#t. (114) 
Solve the Dirac equations by exponentiation, 

Mx) =Vexp[-^ [ X A-(w + ,x-)dw + ]^\x'), <9 + ^ 0) = 0, (115) 

2 J — oo 

if, 2 (x) = Vexp[- 1 -^ [ X A + (x + ,w-)dw-]^\x + ), d~i)f ] = 0. (116) 

2 J— oo 

The point-splitting vector is taken spacelike, e = (O^e 1 ), e^ 1 = ie 1 . Thus, for 
example, ip\ must be replaced by 

ip\(x + e) = $ \x~ -e^Pexpf^ [* + ' A~(w + , x~ - e 1 )dw + }. (117) 

2 J— oo 

The symbols T 3 and stand for path and antipath ordering, respectively. The result 
of the point-splitting differs from the usual one by having integrated (nonlocal) 
terms. The equations for A ± become 

p 2 p 2 rX + Pi A — 

(□ + — )A + -— / ^—{w + ,x-)dw + = 2e^f ) \x-)4>f\x-) (118) 
(□ + —)A~ - — / — (z+, W -)d«;- = 2e^V)V4 0) (* + )- (119) 

Z7T Z7T J-oo CX" 1 " 

The integrated terms here can be removed by taking derivatives with respect to the 
upper limit. Combining the resulting equations leads to 



□ d ■ A = 
25 



(120) 



(□ + -) t^A v = 0; (121) 

71 

thus d ■ A = rj is a massless field and e liV d il A v = £ is a field of mass ej\pK. 
The fields r\ and S (which is the electric field in 1+1 dimensions) are the gauge- 
variant and gauge-invariant degrees of freedom, respectively. Then □ A 1 must 
be a linear combination of d^rj and & v dJL. A 1 must be the convolution of the 
A(x) = — |e(a; )A(x) Green's function with this linear combination plus terms an- 
nihilated by □. The convolution of A (a;) with rj does not exist, because, formally, it 
is / d 2 yA(x — y)r)(y) = J d 2 kexp(—ik ■ x)5(k 2 )fj(k)/k 2 , which is ill-defined. Because 
of this, a new field a that obeys □ a = rj must be introduced. This was first done 
by Capri and Ferrari [f23|j . Thus 

A 1 = Cl d"a + czerdpE + c^rj + c 4 ^ (0) 7 M ^ (0) . (122) 

For the massless case, 

^(°)yy°) = dy, □ = 0, (123) 
where is a free positive-metric scalar field. Using 

[V, v}et = riS, [rj, <x}et = c ari i5, [a, rj\ E T = c av i5, [a, a} ET = c aa i5. (124) 

with the choices c ar} = c aa = and 

/ I 7T / 5"7T 

r = -l, c 1 = T\—, c 2 = ±J — , c 3 = ±W— . (125) 
V 57r v e z V e 

yields a solution in which the vector potential obeys the canonical commutation 

relations. 

The canonical commutation relations in field theory and their predecessors in 
classical mechanics and quantum mechanics are important for many reasons. The 
Poisson brackets in classical mechanics, for example, ensure that the Hamiltonian is 
the generator of time translations. In quantum mechanics, for example, the relation 
[x,p] = ih leads to the uncertainty relation. In quantum field theory, the CCR's lead 
to the free field being a collection of quantized oscillators. In nonrelativistic field 
theories at least, the CCR's imply unitarity[[|. A new feature of the canonical com- 
mutation relations in quantum field theory is that they ensure that the asymptotic 
fields have the proper free commutation relation. (The renormalized canonical com- 
mutation relations will do as well as the original CCR's for this purpose.) For these 
reasons, a solution that obeys the canonical commutation relations is important. 
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7. SUMMARY AND OUTLOOK FOR FUTURE WORK 

The N quantum approach using the Haag expansion has met the test of non- 
gauge theories, including bound states and both spontaneous and dynamical sym- 
metry breaking. Previous work by Amit Raychaudhuri|24| and work presently in 
progress in collaboration with Eli Hawkins ||25|| and with Rashmi Ray and Felix 
Schlumpff^] in both non-gauge theories and in non-confining gauge theories give 
promise of using the N quantum approach to treat bound states in a way that has 
advantages over the Bethe-Salpeter equation. The treatment of confined degrees of 
freedom in quantum chromodynamics remains a goal for the future. 
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